Design and analysis of photonic crystal coupled cavity arrays for quantum simulation 
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We performed an experimental study of coupled optical cavity arrays in a photonic crystal plat- 
form. We find that the coupling between the cavities is significantly larger than the fabrication- 
induced disorder in the cavity frequencies. Satisfying this condition is necessary for using such cavity 
arrays to generate strongly correlated photons, which has potential application to the quantum sim- 
ulation of many-body systems. 



I. INTRODUCTION 

Solving strongly correlated quantum many-body sys- 
tems is a formidable task. One promising approach is to 
mimic such complicated systems using another simpler 
and easily controllable quantum system, as envisioned by 
Richard Feynman [1]. To that end, the first demonstra- 
tion of quantum phase transitions with ultra-cold atoms 
in an optical lattice [2] sparked a significant amount of 
research on quantum simulation with atomic systems [3] . 
Another very promising direction of using photons them- 
selves as the interacting particles has generated consider- 
able interest recently [4] . The main idea of this approach 
is to obtain a correlated "quantum fluid of light" [4] by 
building a coupled network of nonlinear electromagnetic 
cavities. The photons can hop between cavities due to the 
electromagnetic coupling and can repel each other in the 
same cavity due to the intra-cavity nonlinearity. We note 
that such a coupled cavity network exhibits rich physics 
such as topologically protected optical delays [5], even 
without any nonlinearity, although having nonlinear cav- 
ities opens up many more avenues of research. Obviously, 
the optical nonlinearity required for significant repulsion 
at low photon number is very high, and in current tech- 
nology, only 2-level systems (for example, atoms, single 
quantum emitters such as quantum dots (QDs) or super- 
conducting transmon qubits) strongly coupled to a cavity 
provide such strong nonlinearity in the photon blockade 
regime [6-9] . In most of the applications relating to quan- 
tum simulations, one needs to deterministically position 
single quantum emitters in each of the cavities, which 
is very difficult to achieve in the state-of-the-art solid- 
state technology. However, recently several groups have 
demonstrated deterministic positioning of self-assembled 
QDs [10], and the hope is that these site-controlled QDs 
will also perform well within the setting of cavity quan- 
tum electrodynamics (CQED). Another approach would 
be to use a bulk nonlinearity or quantum well nonlin- 
earity, but significantly enhanced by a cavity with high 
quality (Q) factor and low mode volume [11, 12]. We 
note that such a platform consisting of coupled nonlin- 
ear cavities is useful not just for the quantum simulation, 
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but also for quantum error correction [13] as well as for 
classical optical signal processing [14] . 

Although plenty of theoretical proposals for simulating 
interesting physics in such a coupled cavity array (CCA) 
are present in the literature, the experimental progress 
in that direction is rather limited. As one needs to have 
many cavities for this operation, a solid-state system is 
obviously an ideal choice. However, due to imperfect 
nano-fabrication solid-state cavities have inherent dis- 
order, resulting in different resonance frequencies than 
the cavities were originally designed for. Such disorder 
might limit the utility of CCAs for quantum simulation. 
However, in a recent paper it is argued that as long as 
the coupling strengths are much larger than the disor- 
der, the CCAs can be used for quantum simulation, and 
it is shown that microwave transmission line cavities for 
circuit QED satisfy this condition [15]. 

In this paper, we demonstrate high-Q 2-D CCAs based 
on photonic crystals fabricated in GaAs with embedded 
high denisty self-assembled InAs QDs. Although a pair 
of coupled cavities, also known as a photonic molecule, is 
well studied in the literature [16-18], relatively little lit- 
erature exists for CCAs. A 2-D CCA of photonic crystals 
in GaAs (with multiple quantum wells as active materi- 
als) has been studied previously for increasing the out- 
put light intensity from nano-lasers or slowing down light 
[19-22], but the Q-factors of the cavities were too low to 
identify individual cavity modes. A long chain of high-Q 
coupled cavities has been studied in silicon [23], but the 
physical phenomena observable in such a 1-D chain are 
rather limited. While a 1-D array [24] has been studied 
as a platform to simulate the physics of Bose glass [25] , 
and Tonks- Girardeau gas [26], a 2-D array is a more suit- 
able candidate for simulating many other systems includ- 
ing topologically non-trivial states such as the fractional 
quantum Hall state [27-29]. 



II. SPECTRA OF COUPLED CAVITIES 

In our experiment, we employ an array of linear three- 
hole (L3) defect photonic crystal cavities, typically stud- 
ied in single QD-cavity QED experiments [30, 31]. The 
fundamental mode of such a cavity is linearly polarized 
in the direction orthogonal to cavity axis; in our pro- 
posed CCA geometry all the cavities have parallel axes, 
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FIG. 1: (color online) Scanning electron micrograph (SEM) 
images of CCA with (a) 4 cavities, (b) 9 cavities and (c) 16 
cavities. The simulated electric field profiles for each of the 
two super-modes of the coupled cavities are shown: (d),(e) 
for 60° coupled cavities; (f),(g) for laterally coupled cavities; 
(h),(i) for vertically coupled cavities. 
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FIG. 2: (color online) The PL spectra of the CCA for (a) 
4; (b) 9 and (c) 16 cavities. We can clearly identify all the 
cavity array modes. We focus on several specific separations 
between the CCA modes labeled Ai through A6 in the plots 
and perform statistical analysis. We also observed several low- 
Q modes at long wavelengths for several cavity arrays, as can 
be seen in part (a). These modes are not the actual coupled 
cavity modes under study, which is confirmed by monitoring 
the resonance frequencies of single (uncoupled) L3 cavities 
fabricated in the same chip. 



and their modes have the same polarization. The pho- 
tonic crystal CCA is fabricated in a 164 nm thick GaAs 
membrane (with self-assembled InAs QDs embedded at 
a depth of 82 nm from the surface) using electron-beam 
lithography and reactive ion etching [30]. Scanning elec- 
tron micrographs of the fabricated structures are shown 
in Figs. la,b,c. Three different CCAs are designed con- 
sisting of 4, 9 and 16 cavities. These cavities are coupled 
to each other by three different coupling strengths de- 
pending on the relative orientation and separation of two 
cavities. When two cavities are coupled at an angle of 60° 
(Figs. ld,e) the coupling strength t is strongest; for ver- 
tically stacked coupled cavities (Figs. lh,i) the coupling 
strength J\ is smaller than t] and for horizontal coupled 
cavities (Figs. lf,g) the coupling strength J2 is much 
smaller than t and J\ (the difference in coupling strengths 
is a result of the different radiation patterns of the cavity 
modes, and their different overlaps in various directions). 
From the finite difference time domain (FDTD) simula- 
tions we can calculate the field profiles of the coupled 
cavities (Figs. 1 d-i) and estimate the coupling strengths 
from the separation of the super-modes in the simulated 
spectra, assuming cavity operation in the range of QD 
emission (~ 900 — 930 nm). For a hole radius r varying 
from 70 nm down to 50 nm, with photonic crystal lattice 
constant a = 264 nm, we find that £/27r ~ 0.8 — 1.3 THz; 
Ji/2tt - 0.4 - 0.8 THz and J 2 « t,J^ 

We characterize the resonances of the coupled cavity 



array by photoluminescence (PL) studies, where the large 
density of embedded QDs acts as an internal light source. 
Fig. 2 shows the PL spectra obtained from the CCAs. 
The spectrum for different CCAs are taken from different 
e-beam doses, and hence correspond to slightly different 
hole radii. The quality factors of the observed modes 
are ~ 1000 — 3000, and all the modes are linearly po- 
larized with similar polarization axis. We note that the 
set of higher Q-factor resonances in Fig. 2 are identi- 
fied as the coupled fundamental modes of the L3 cavi- 
ties, shown in Fig. 1. These modes are not necessarily 
in the same wavelength range for different sizes of the 
arrays, as the structures were defined during the fabrica- 
tion process with different doses in e-beam lithography, 
and thus photonic crystals have different parameters. We 
also point out that the number of modes observed in PL 
should be the same as the number of cavities in the CCA, 
irrespective of whether the cavities are coupled or not 
(assuming any degeneracy is lifted due to fabrication im- 
perfection). Without coupling between the cavities, the 
observed modes would be randomly placed and no spe- 
cific order between the modes should be observed. On 
the other hand, in the presence of the coupling between 
the cavities, the cavity modes are expected to be spaced 
at a specific order determined by the coupling strengths. 
However, due to the disorder introduced during the nano- 
fabrication process, the exact distribution of the cavity 
resonance frequencies will be perturbed. Hence from a 
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statistical study of the differences in the cavity resonance 
frequencies we can estimate the ratio between the cavity 
coupling strengths and the disorder in the cavity reso- 
nances. We note that one could instead estimate the 
disorder in the cavity resonances from the actual cavity 
frequencies, and not the differences. However, cavities 
written on different parts of the chip are more suscepti- 
ble to fabrication variation, and might suffer an overall 
frequency shift. Thus, the mode separations provide a 
better measure of the disorder present within each CCA, 
while allowing us to gather statistics from several CCAs 
for comparison. 

We find a consistent order between the modes of dif- 
ferent CCAs (Fig. 2), indicating the cavities are coupled. 
Next we analyze all the separations between the subse- 
quent cavity modes. In order to do this, we fabricated 
~ 30 copies of each of the three types of cavity arrays, 
and calculated the mean /i and standard deviation er of 
all these mode separations. We note that the mean \i of 
the mode separations is a combination of the coupling 
strength and the disorder, whereas the standard devia- 
tion er of the mode separations depends mostly on the 
disorder, as explained earlier. To elaborate further, we 
can consider the simple example of a photonic molecule 
(two coupled cavities), where the observed separation A 
between two modes is \/Aq + 4 J 2 with Ao being the ran- 
dom bare detuning between the cavities due to fabrica- 
tion imperfection and J being the the coupling strength 
[16]. We assume that the bare detuning follows a Gaus- 
sian probability distribution with zero mean and stan- 
dard deviation <tj, i.e., the probability of having a de- 
tuning Ao is 

A 2 

Pr(Ao) = ^=^e"4 (1) 

Under such a Gaussian approximation [15], we find that 
the mean \i of the mode separation A (we consider 

the absolute value of the separation) is \i = \J~^erf if 

there is no coupling (J = or (i.e, erf /J » 1) and 
/i = 2J-\-a < j/4J-\-0(aj) if the disorder is weak compared 
to the coupling i.e., erf /J << 1). The standard devia- 
tion er for the mode separations is er ~ Of without any 
coupling (J = 0) and er ~ 0(a'j/J) when erf /J « 1. 
Similar analysis can be performed for CCAs with more 
than two cavities, although the expressions for the mean 
and standard deviation become complicated, and a sim- 
ple closed form expression is difficult to obtain. Never- 
theless, as seen for the photonic molecule, the ratio of 
the standard deviation to the mean gives us an estimate 
of the relative contribution of the disorder and the cou- 
pling to the mode separations. In our fabricated CCAs, 
we find that the ratio a/fi « 1 for almost all the mode 
separations, indicating the presence of strong inter-cavity 
coupling; o therwise , for no coupling, the ratio would be 
equal to \/tt/2 — 1 ~ 1. Table I shows the data for spe- 
cific separations (Ai — >• Aq) between the cavity modes 
in the cavity array spectra in detail. We note that all 



TABLE I: The mean mode separations (/x), and standard de- 
viation (a) measured over ~ 30 cavity arrays, with similar 
hole radii (see Fig. 2 for definition of the separations). 
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FIG. 3: (color online) Numerically-calculated eigen-spectra 
of the coupled cavities: the eigen- frequencies as a function of 
the disorder standard deviation af for (a) 4, (b) 9 and (c) 
16 cavities in the arrays. The spacings between the cavities 
are the same as the structures shown in the SEM images in 
Figs. 1 a-c. The difference in the subsequent eigen-values are 
shown as a function of 07 for (d) 4, (e) 9 and (f) 16 cavities 
in the arrays. We note that the mode separations increase 
linearly with increasing a/, when a/ is much greater than 
the coupling strengths, as found in the theory from a simple 
photonic molecule. Insets magnify the region of low disorder 
and we identify the mode separations Ai — >> A6. 



the separations are not equally influenced by the cou- 
pling strengths as seen from the numerical simulations 
presented below (Fig. 3), and the chosen separations 
(indicated in Fig. 2) are the ones that are most heavily 
influenced by the coupling strengths. 
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III. ESTIMATION OF COUPLING AND 
DISORDER 

Using the coupling strengths derived from FDTD sim- 
ulations (*/2tt = 1.2 THz, Ji/2tt = 0.8 THz, J 2 « 0), we 
calculate the eigen-states of the CCA by diagonalizing 
the Hamiltonian 'H: 

H ^Aiajfli + ^^(aja^- + ata*) (2) 

where, Ai is the resonance frequency of the i th cav- 
ity due to fabrication imperfection, and g^j is the cou- 
pling strength between the i th and j th cavity. The cav- 
ity frequencies A^'s are randomly chosen from a Gaus- 
sian distribution with zero mean and standard devia- 
tion (jf (given by Eq. 1). We assume that the coupling 
strengths are constant and do not depend on the disor- 
ders. The mean values of the eigen-frequencies (averaged 
over ~ 10000 instances) are plotted in Figs. 3a,b,c as a 
function of increasing cry. We observe that the relative 
detunings between the modes follow a specific order when 
the disorder is small. However, with increasing disorder 
any specific order between the cavity modes disappears. 
This can be observed more clearly in Figs. 3d,e,f, where 
the differences in the mode frequencies are plotted as a 
function of (jf. We note that the differences become sim- 
ilar, and increase linearly with <jj, as found in the simple 
photonic molecule model previously. We also note that 
several modes are spaced very closely at weak disorder, 
indicating a lesser contribution of the coupling to such 
mode separations (inset of Figs. 3d,e,f). On the other 
hand several detunings between the modes are large com- 
pared to others (denoted by Ai — » A 6 ) signifying a large 
contribution from the coupling strengths to the mode 
separations. We observe that the relative positions of 
the cavity modes match qualitatively with our experi- 
mental findings. We can identify the specific separations 
Ai —> Aq between the modes (Fig. 3). Clearly, the fab- 
ricated structures are in the regime where the coupling 
strengths are greater than the disorder. This regime is 
magnified in the inset of Figs. 3d,e,f; and the mode sep- 
arations Ai — » A6 are identified. 

Finally, as a further proof of the fact that the detun- 
ings between the observed cavity array modes are mostly 
due to the coupling between the cavities, and not the 
disorder, we repeated the fabrication of sets of ~ 30 cav- 
ities for different values of air-hole radius for all three 
types of CCAs. A decreasing trend in the separation is 
observed with increasing hole radius (Fig. 4). A similar 
trend is observed in simulation for a photonic molecule 
with diagonally placed cavities, as a function of the hole 
radius (inset of Fig. 4). Such a trend also indicates that 
the separations are mostly due to the coupling between 
cavities, as a detuning due solely to disorder would have 
a much weaker dependence on the photonic crystal hole 



size. The decrease in the mode separation with increase 
in the hole radius can be explained by the increase in 
the photonic band gap size with increasing hole radius 
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FIG. 4: (color online) The mode separations Ai — > A6 as 
a function of the photonic crystal hole radius. A decreasing 
trend in the separation is observed with the increasing hole 
radius (the photonic crystal lattice periodicity a is 264 nm). 
For comparison, the inset shows the numerically simulated 
(FDTD) coupling strength between two cavities placed diag- 
onally as a function of the hole radius. 

(and thus larger reflectivity of the mirror layers separat- 
ing cavities, which reduces the cavity couplings). 

IV. CONCLUSION 

We show the signature of large coupling strengths be- 
tween photonic crystal cavities, in a coupled cavity ar- 
ray fabricated in GaAs containing InAs QDs. We ob- 
serve that the coupling strengths are significantly larger 
than the disorder introduced during the nano- fabrication. 
Satisfying this condition is necessary for employing such 
cavity arrays in quantum simulation with correlated pho- 
tons, although the challenge of achieving a nonlinear ity 
in each cavity still remains open. 
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